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ABSTRACT 

Aims. The appearance of folding ion rays in cometary comae is still not very well understood, so our aim is to gain more insight into 
the role of the local solar wind in the formation of these structures. 

Methods. Comet C/2004 Q2 (Machholz) was intensively monitored during its closest approach to Earth (January 2005) with the 
CCD camera Merope mounted on the Flemish 1.2m Mercator telescope, in three different bands (Geneva U and B and Cousins I). 
Spectacular ion rays, thin ionic structures rapidly folding tailward, were recorded in the U band during one night, January 12th. 
Results. Data from the SOHO satellite that was extrapolated corotationally to the position of the comet showed that the ion rays were 
formed during a sudden change in the in-situ solar wind state. We were able to succesfuUy correlate a high-speed solar wind stream 
with the appearance of folding ion rays. 

Conclusions. To our knowledge, this is the first clear observational evidence that folding ion rays in cometary comae are produced 
by a sudden change in the local solar wind state. 
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1. Introduction 

A comet can be considered to consist of a small nucleus 
(~10km) surrounded by a dusty gas cloud, dubbed the coma 
(~10^km). On even larger scales, the coma is embedded in the 
halo, extending up to ~10^km. Two other typical components 
of comets in the vicinity of the Sun are the dust tail, consist- 
ing of relatively large neutral grains, and the ion tail, made of 
ionized gas. Basically, the two tails are seperated because the 
outward solar radiation pressure is higher than the inward grav- 
itational pull. |Biermann| ( |1951| l showed that radiation pressure 
alone could not be responsible for the outward force on the ions. 
In this way, the morphology of comets was one of the first proofs 
for the existence of the solar wind. 

Thus, variations in the solar wind result in a variable ion tail. 
Conversely, through certain properties of the cometary plasma 
structures, it may be possible to deduce valuable information on 
the solar wind. For example, a recent survey of Chandra comet 
observations has demonstrated that X-rays provide a quantita- 
tive diagnostics of the interacting solar wind (Bodewits et al.| 



2007[ l. Other properties that have been exploited are tail dis 



connection events (DEs); see for example Niedner et al. ( 1978 1, 
|Snow et"aL] ( |2004[ ), and the recent tail rip-ofF from comet 2/P 
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Encke ( |Vourlidas et al.|2007] l and statistical studies of the aberra- 
tion angle of the main tail axis ( Brandt et al.|1973j ). This sharply 
contrasts with the study of cometary folding rays, as they ap- 
pear closer to the nucleus and thus require a higher spatial res- 
olution. Current theories suggest that folding rays are tailward 
extensions of ion pile-up regions (receding envelopes) on the 
sunward side of the nucleus, but it is still not well known what 
exactly causes this phenomenon. Possible explanations include 
external origins: some studies invoke the crossing of the comet 
with IMF (interplanetary magnetic field) sector boundaries, oth- 
ers denote the IMF as just a source mechanism and put more 
emphasis on the cometary ionosphere. Still other studies claim 
they have a completely internal origin and are dependent on the 
nucleus' ion production rate (see Ip 2004' for a review). This let- 
ter supports the idea that the solar wind is largely responsible for 
the appearance for cometary folding rays, as we found a clear 
indication that a fast cold wind (originating from coronal holes) 
triggers tail rays. 

We observed comet C/2004 Q2 (Machholz) while it encoun- 
tered a change in the solar wind regime. Such observations are 
rare, as there are a lot of variables that have to be set: firstly, 
we must have accurate solar wind data and detailed images 
of the comet with an adequate spatial and temporal resolution. 
Secondly, the observed comet has to be relatively close to the 
Earth, preferably in or near the ecliptic, to be able to extrapolate 
solar wind information to the position of the comet (as the so- 
lar wind is highly variable in both time and place). Furthermore, 
a solar event has to occur during the observing run and the ef- 




Fig. 1. Azimuthal renormalised frames of 6 of 8 observations in the U band of comet C/2004 Q2 (Machholz) on 2005 Jan. 12, 
together with three observations of the night before (upper left), and one on the night after (lower right). The trails on the images are 
due to background stars and indicate the comet's motion. All images taken on the 12th show clear ion rays, while on the 11th, only 
the typical sharp ion tail is visible; no ion rays appear. In the right lowerhand corner, it is clear that the ion rays gradually disappear. 
All images have the same scale and color coding, the direction of the Sun is upwards, and the field of view radius is approximately 
^SxlO'^km. 



fects have to be detected. It is quite a challenge to prepare for 
such an observation run, as many of these variables are not pre- 
dictable on a long term and thus constant monitoring of solar 
wind data is necessary. In this light, the passage of comet C/2004 
Q2 (Machholz) in January 2005 was an exceptional event since 
all of these conditions were set. 



2. Observations, reduction, and image 
enhancement 

We obtained about 181 CCD-observations of comet C/2004 Q2 
(Machholz), taken in three different broadband filters, namely 
Geneva U (61 observations of 600 sec integration), Geneva B 
(61 observations of 60 sec), and Cousins I (59 observations 
of 60 sec). The images were taken with the Merope camera 
mounted on the 1.2m Mercator telescope on La Palma, Spain, 
between the 7th and 15th of January 2005. The camera consists 
of a 2158 X 2044 pixel CCD with a field of view of 6.'5 x 6.'5 
and a resolution of 0."19pix"'. This translates to ~50kmpix"' 
at the position of the comet. The reduction of the CCD images is 
described in detail in Reyniers et al. ( 2008 1. 

The overwhelming brightness of the coma makes it difficult 
to study relatively faint structures such as ion rays. For a spher- 
ically symmetric, optically thin coma with a stationary outflow 
of dust grains, the projected brightness profile in the plane of 
the sky equals a 1 /r distribu tion, with r the projected radius in 
the sky (Wallace et al. 1958). The inner coma (r< lO'^km) is 
thus very bright in comparison with other features. The situation 
changes beyond the outer parts of the coma, where the overall 
density of the gas and dust becomes too small to reflect much 
light, and the accumulation of ions in certain directions now be- 
comes relatively high. Even in the inner coma, these features are 
present in the signal, but the signal has to be enhanced to make 
them visible to the naked eye. Many techniques have been de- 
veloped (see e.g. [Schleicher & Farnham|2004[ for an overview). 



but each technique suppresses certain features while bringing 
out others, so they have been carefully examined and tested. In 
general, these techniques include dividing and subtracting the 
frames by empirical or theoretical (sunward) surface brightness 
profiles. The downside here is that we lose information on az- 
imuthal structures like shells (Schulz 1991 ). 

Other techniques include subtracting/dividing with ro- 
tated/translated copies, azimuthal renormalisation around the 
nucleus, and filtering with diff'erent kernel sizes, such as the 
Larson-Slaughter filter ( Larson & Slaughter|1992[ l and Larson- 
Sekanina filter ( |Sekanina & Larson|1984| l. We chose to subtract 
an azimuthally renormalised profile via 



/(r,0) = /o(r,0)- 



z- 



/(r) 
{r) 



where n{r) is the number of pixels on a distance r from the nu- 
cleus. To avoid artifacts due to the presence of field stars, it is 
usually better to replace the average by a median. In our case, 
the technique of azimuthal renormalisation has been shown to 
be both the fastest and qualitatively best procedure to enhance 
highly radial features such as ion rays, because it produces vir- 
tually no artifacts. 



3. Folding ion rays 

If we take a look at our dataset, we see that the observations in 
the U-filter display a very variable behaviour in the tailward di- 
rection, while the images in the I-band do not seem to change 
at all. This is a clear indication that the I-filter focuses on the 
cometary dust, while the U-filter shows ionic structures. The U- 
band is dominated by fluorescent emission from mo lecules like 



CO+ (3787 A), N+ (3914 A), and OH+ (3600 A) (Lutz et al 



|1993). For example, on 11 January, the typical sharp ion tail 
of the comet is clearly visible. The most remarkable event is, 
however, the appearance of thin ion rays on 12 January. On that 
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Fig. 2. Polar representation of images taken 2005 Jan. 12 at 
22:50 and 23:06 UT. The comet's nucleus is to the right. In the 
upper section of both images, two ion rays are visible (indicated 
by tickmarks). These rays are folding along the tail's axis with 
an angular velocity of ~ 0.9 degrees per minute. 



date, the observations differ significantly compared to the images 
taken on previous and consecutive days; before they are not vis- 
ible, and afterwards, they seem to disappear gradually (Fig.[T]i. 
Further investigations of the frames reveal a very symmetric pat- 
tern for about all observations on 12 January. The ion rays are 
symmetric around the tail axis, although possibly small devia- 
tions cannot be excluded. 

On images taken at 22:50 UT and 23:06 UT (16 minute in- 
terval), the evolution of two ion rays can be traced. It is visually 
clear from their polar representation (Fig. |2]i the rays are fold- 
ing along the tail axis of the comet. To derive an estimate of the 
angular velocity with which the rays converge, we integrated the 
flux radially over an interval of 70 pixels, starting at a distance of 
530 pixels from the optocentre. This interval was chosen because 
of the high signal and the absence of stars. This way the signal 
was significantly enhanced, which made it possible to represent 
the position of the ion rays by the local maxima of a polynomial 
fit through the integrated flux as a function of the azimuthal co- 
ordinate. We concluded that these maxima folded around the tail 
axis with an angular velocity of ~ 0.9 degrees per minute. To es- 
timate the order of magnitude of the physical velocity associated 
to this value, we calculate that, at a distance where 1 pixel spans 
1 angular degree (which is at about 3000 km from the nucleus), 
the velocity of the structure is ~235kms"'. This is too high to 
be of cometary origin, whereas it has the same magnitude as the 
solar wind speed. This brings us to a closer examination of the 
solar wind conditions at the position of the comet. 

4. Connection with solar wind through SOHO and 
ACE data 

SOHO and ACE are located at the first Lagrangian point LI be- 
tween the Sun and the Earth (ca. 0.01 AU from Earth). One of 
the many purposes of the SOHO satellite is to measure the so- 
lar wind velocity, which is done by the Proton Monitor on the 
CELIAS instrument ( [Hovestadt et al.||1995| ). Together with the 



measurements of the solar wind proton density taken by the same 
instrument, this gives us the solar wind flux as well. 

From a first look at the SOHO data (Fig. [3] upper panel), 
a clear discontinuity is seen in the solar wind velocity on 12 
January at LI. We cannot immediately presume that this event 
also happened at the comet C/2004 Q2 (Machholz) position, be- 
cause the solar wind has a finite velocity that is variable in both 
space and time as well. First of all there are two regimes of 
solar wind, namely the "slow" and "fast" regimes, with veloc- 
ities around 400 km s"' and 700 km s"', respectively. Also, it is 
known that the solar wind can be very different at different lat- 
itudes. We can safely assume that the sudden rise in the wind 
velocity is due to a corotating interaction region (CIR) originat- 
ing from a coronal hole, because there was a peak in the proton 
flux just before the wind velocity rose and stayed up for a day or 
two (Fig.|3] middle panel), which is typical of a CIR (Niednerj 

We used the measurements of the O /O -ratio from 
ACE/SWICS dGloeckler et al.|1998| l to detemiine the nature of 
the solar wind streams (i.e. originating from coronal holes, coro- 
nal mass ejections (CME's), or simply slow solar wind). Oxygen 
was chosen because it is the most abundant heavy ion. The ra- 
tio of heavy highly-ionized ions O^^/O^^ gives us an idea of 
the temperature of the wind, therefore also of its origin. The 
low observed 0''^/0''^ ratio implies a low-freeze in tempera- 
ture, suggesting a cool origin on the sun. A ratio around 1 % is 
indeed in good agreement with values typical of a coronal driven 
wind (Zurbuchen & Richardson 2006) (Fig. [3]l. Moreover, the 
SOHO/LASCO instrument detected no CME's heading for Eaith 
at that time (http://lasco-www.nrl .navy.mil/halocme.html). 

The velocity of a CIR has both a radial as a corotational com- 
ponent due to the rotation of the Sun. Thus, if we want to extrap- 
olate the SOHO/CELIAS measurements to the position of comet 
C/2004 Q2 (Machholz) during the observation run, we have to 
make certain that the comet is not situated too far away from LI, 
that it is not too high above or below the ecliptic, and that we 
take both components of the outflow velocity into consideration. 
The first two conditions are immediately set (Fig. |4|i. Next, the 



actual extrapolation is done by using the formula (Neugebauer 
let al.|2000| ) 

Af = Mrad + Af,„, . (1) 

The first term on the right hand side of equation ([T]) equals 
Ar/v0 with Vq the velocity of the solar wind and Ar the dif- 
ference between the heliocentric distances of LI and comet 
C/2004 Q2 (Machholz). We can set the velocity of the comet 
to vj. = because <K Vq. The second term on the right hand 
side can be rewritten as Af™, = A/o„/Qq, with A;„„ the longi- 
tudonal differences between LI and the comet's position, and 
Q.Q - 2.9 X 10 ''rads ' the solar rotation frequency. For exam- 
ple, when the ion rays become visible for the first time, we cal- 
culate A? ~6 h. We emphasise that this method does not provide 
us with exact values of solar wind parameters at the comet's lo- 
cation, as it does not account for any magnetohydrodynamical 
phenomena, such as shocks and CMEs. That is why the extrap- 
olation of the SOHO data does not result in a smooth function, 
but sometimes gives two predictions of the solar wind velocity 
at the same moment: wind particles travelling faster but origi- 
nating later than slower particles may, according to our approxi- 
mations, aiTive at the comet's position at exactly the same time, 
rather than result in a shock front. 

Taking these limitations into consideration, we can conclude 
that the change from a slow to fast solar wind regime, in the 
form of a corotating interaction region detected on 12 January in 
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Fig. 3. Upper panel: Measurements of the solar wind veloc- 
ity from SOHO/CELIAS (black line) and a corotational map- 
ping (red line) at the location of C/2004 Q2 (Machholz). Middle 
panel: The proton flux from SOHO/CELIAS Proton Monitor at 
LI. The solar event on 12 January was due to a CIR, which 
is typically preceded by a sudden peak in the proton flux. 
Lower panel: Ratio of 6'^ jO^^ from ACE/SWICS. The red 
line represents the results from an empirically deduced formula 
( |Zurbuchen & Richardson||2006| l, for which a low value (hor- 
izontal red line) indicates a fast, cold solar wind. The vertical 
black line denotes the first observation on 12 January. 
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Fig. 4. Configuration of C/2004 Q2 (Machholz), the Earth, 
and Sun on 2005 Jan. 12, when the solar wind regime 
changes suddenly from slow to fast, and ion folding rays 
are visible. All orbits are seen from above and projected 
onto the ecliptic. Note that Machholz' orbit has an incli- 
nation of ~39deg with respect to the ecliptic (data from 
|http : //ssd ■ jpl ■ nasa. gov/horizons . cgi| l. 



LI, is expected to have arrived at C/2004 Q2 (Machholz) on the 
same date the ion rays are visible in our images. 

5. Discussion 

Overall, there are three possible explanations for the observed 
features. First, they could be caused by eruptions from the nu- 



cleus. This is highly unlikely, because the ion rays are too numer- 
ous, too symmetric, always located in the tailward hemisphere, 
and they change their position too rapidly (e.g. the images taken 
on 22:50 UT and 23:06 UT, Fig.[T). Second, the tail rays could 
look like elongations of the shock front. This is not very li kely ei- 
ther: considering a gas production rate Q of 2.6x 10^'s"' (Bonev 
|et al.|200 6), the shock front is expected to be situated at a dis - 
tance Rs > 90 x 1 0^ km from the nucleus ( Wegmann et al. 2004 1, 
which is just outside the field of view, but the features are likely 
to be within the bow shock. A second counterargument is that a 
shock front is only visible in X-rays ( Wegmann et al.|[l998 l or 
by in-situ measurements (Fuseli er et al.|1991[ ). A third explana- 
tion invokes discontinuities in the solar wind (jlp 2004). This is 
the most likely explanation since such a discontinuity was de- 
tected by SOHO around the time of the observation. Apparently, 
a discontinuity in the solar wind, caused by a CIR, is a good can- 
didate for triggering the appearance of cometary folding rays. 



6. Conclusion 

Comet C/2004 Q2 (Machholz) reached a visual brightness of 
~3.5 mag. at its closest approach to Earth on 2005 Jan. 5, mak- 
ing it one of the most spectacular astronomical events of that 
year. The comet was visible during a significant part of the night, 
making it ideally suited to performing a dedicated ground-based 
monitoring. The first objective of our intensive monitoring was 
to derive a rotation period ( [Reyniers et al.|2008) l. In addition to 
this, coma structures were also studied in three different bands. 
During one night, 12 January, clear ion rays were detected in 
the U band images after enhancement. Due to the comet's po- 
sition relatively close to the Lagrangian point LI and close to 
the ecliptic, we were able to extrapolate the solar wind veloc- 
ity data from SOHO to the comet's position. It turned out that 
the development of the ion rays coincide perfectly with a sudden 
change in the local solar wind from the "slow" (~ 400 km s"') to 
the "fast" (~ 700km s"') regime. To our knowledge, this is the 
first time that the relatively rare event of folding ion rays in the 
coma of a comet could be correlated with a change of the in-situ 
solar wind state. 
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